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Motivated by the transmission of topological surface states through atomic scale steps, we study 
the transport of gap less Dirac fermions on hyperbolic surfaces. We confirm that, independent of the 
curvature of the hyperbolae and the sharpness of the corners, no backward scattering takes place 
and transmission of the topological surface states is completely independent of the geometrical 
shape (within the hyperbolic model) of the surface. The density of states of the electrons, however, 
shows a dip at concave step edges which can be measured by an STM tip. We also show that the 
tunneling conductance measured by a polarized scanning tunneling probe exhibits an unconventional 
dependence on the polar and azimuthal angles of the magnetization of the tip as a function of the 
curvature of the surface and the sharpness of the edge. 
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I. INTRODUCTION 

Interest in topological insulators (TIs) continues to re- 
main high ever since their prediction pQ and discovery [2]. 
The hallmark of these materials is their topologically 
non-trivial nature, due to which they have gapless lin- 
early dispersing edge states, although they are insulating 
in the bulk [3]. Such edge states are familiar in the two- 
dimensional quantum Hall system [4 , but in TIs, they 
arise even in higher dimensional systems and even in the 
absence of time-reversal symmetry breaking. 

For two-dimensional TIs, the edge states consist of two 
counter-propagating modes, with opposite spin projec- 
tions. The correlation of the spin of the electron with its 
direction of motion is the central feature of these states. 
The existence of these states have been verified exper- 
imentally by measurements in the HgTe quantum well 
structures [2]. There have also been several theoretical 
studies [SHJ], studying consequences of the spin projec- 
tion and the helical nature of the edge states, although 
there has been no direct measurement of the spin of the 
edge states. 

The three dimensional TIs come in two classes, strong 
and weak, with even and odd number of Dirac cones on 
its surface. For strong TIs, the existence of an odd num- 
ber of Dirac cones is topologically protected and it has 
been found that the surface states of strong TIs such as 
HgTe, Bi2Se3 and BizTe^ can be described by a single 
Dirac electron. These surface states are robust against 
perturbations that do not break time-reversal symme- 
try. The electrons also have the feature that their spin is 
'locked' to the momentum of the electron, which is what 
leads to a complete absence of back-scattering, because 
the spin of the electrons at momentum k is orthogonal 
to the spin of the electrons at momentum — k. The com- 
bination of topological stability and absence of backward 
scattering leads to the prediction that the surface states 
of TIs wrap the surface of the TI and are impervious 
to the existence of surface defects. This has been con- 
firmed by experiments which not only have confirmed the 
absence of backward scattering [8j [9], but more recently 



have also shown full transmission through atomic step 
edges[10 . This would imply that transport through these 
surface states is independent of the geometrical shape of 
the surface. 

With a few notable exceptions [TTHT3] , most of the the- 
oretical work on TIs has been restricted to planar sur- 
faces and straight line edges. For instance, the original 
derivation of the edge states was for samples with a sin- 
gle straight edge. This was extended to the case for finite 
strips where it was shown that there was some interfer- 
ence between the two edges [14 . However, there had been 
no direct study of what happens when two sharp edges 
meet one another. Since spin projection is tied to the 
direction of motion, it is also not clear how the spin cur- 
rent changes if there are sharp edges. Similarly, although 
there have been some studies [T5HT8] of what happens at 
junctions of two 3DTIs, the extension of the gapless Dirac 
state over curved surfaces has yet to be demonstrated. 

A step in this direction was recently taken by Takane 
and Imura[T9]. who introduced a hyperbolic system to 
treat the 90° step edge, and by introducing appropriate 
curvilinear coordinates, they could show that no reflec- 
tion takes place at the 90° step edge, and transmission 
was perfect, although there was a sharp change in the ex- 
pectation value of spin in the close vicinity of the step. In 
this paper, we generalise their work to step edges of arbi- 
trary angles, and show that independent of the curvature 
of the hyperbolae and the sharpness of the corners, no 
backward scattering takes place and the transmission of 
the topological surface states is completely independent 
of the geometrical shape (within the hyperbolic model) 
of the surface. Moreover, we study how the density of 
states (DOS) and the spin DOS behave as a function 
of the curvature and the sharpness of the edge of any 
sample. We find that the DOS shows a dip at the con- 
cave edges of the sample. We also compute the tunneling 
conductance measured by a polarized scanning tunneling 
microscope, as a function of the curvature of the sur- 
face and the sharpness of the edges and show that the 
STM conductance has a non-trivial dependence on the 
curvature angle (j) and an unconventional dependence on 
the polar and azimuthal angles of the tip, which are not 
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displayed by a planar TI surface. 

II. THE MODEL AND ANALYSIS 

We start with the continuum model of a strong 
anisotropic TI in 3D, given by the Hamiltonian [20 
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The parameters A x1 A yi A z and m(p) can be determined 
for specific materials by comparing with the ab initio 
calculations of the effective model of 3DTI[20] . Here, 
the 2x2 matrix represents the spin degrees of free- 
dom and the orbital degrees of freedom are represented 
by the Pauli matrices 7$. The mass term m(p) = 
mo + m2(p x + P 2 y + P 2 Z ) in this 3D Dirac Hamiltonian 
is momentum dependent. For the topological insulator, 
we need to take mo > and rri2 < 0. Assuming the 
system to be translationally invariant in the x direction, 
so that p x is a good quantum number, it is straightfor- 
ward to derive the surface states on either the xy plane 
or the xz plane. However, it is not obvious what hap- 
pens close to the corners. Takane and Imura[T9] studied 
the question of whether or not there exists reflection at 
the corners for non-zero p x by assuming that the surface 
could be represented by a rectangular hyperbolic model 
using curvilinear coordinates. In this paper, we assume 
that the surface of the TI can be an arbitrary hyperbolic 
surface as shown in Fig. (la). Note that we are repre- 
senting a concave surface - the TI is in the shaded region 
as shown in the figure. The curve in Fig. (lb) can be 
described by the equation 
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where the angle between the two asymptotic surfaces is 
governed by the (curvature) parameter a and the sharp- 
ness at the edge is governed by the (sharpness) param- 
eter R. The 3D TI is translationally invariant in the 
x-direction. A separable coordinate system can now be 
defined as follows - 



y = -rcos</> + ?/o(<^ 
z = -r sin0 + zo (</>), 



(3) 



where yo and zo define the curve and are given by 
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We draw normals to the surface as shown in Fig. (lb) 
and define <p as the angle between the normal and the y- 
axis, and define r as the distance (in the direction shown 
by the arrow) from the surface. It is easy to see that (j) 
ranges from tan -1 (a) to 7r — tan -1 (a) and r ranges from 




FIG. 1. The surface is shown (above) in Fig. (la) with the TI 
filling the concave part and the curvilinear coordinate system 
is shown (below) in Fig. (lb) as described in the text. 



— oo to oo. The Jacobian of the transformation can be 
easily found to be given by r + f ((/>), where 
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We can now define the mass term m(p) in terms of 
p r ,P(f) and p x - i.e., m(p) = m r + m^ + m x with 

m r = m - m 2 [<9 2 + (r + /) _1 S r ], 
m = -m 2 [(r + /)- 2 (S 2 - (r + /)-(^/)5 ), 
and m x = rri2p 2 : . (6) 

The bulk Hamiltonian H = m(p)r z -\-(A x p x o- x -\-Ayp y (jy-\- 
A z p z a z )r x can now be written in terms of the curvilinear 
coordinates (r, 0, x) as 
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and By = 



We have used B z = A z (r + /) _1 cos 
A y (r + /) _1 sin((/>)<90 in above. Note that for a rect- 
angular hyperbolic surface with a = 1, our surface and 
expressions are equivalent to to those given in Ref.[19 , 
albeit rotated by 7r/4 about the x-axis. 

To derive the effective 2D surface ((/>, x) Hamiltonian 
from the above bulk Hamiltonian, we need to first solve 
the radial equation H r \ip) = E r \ip). Following Refs. [TT| 
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[19], we obtain a solution of the form = e~ Kr \u), where 
ft -1 measures penetration into the bulk, provided that 
we assume that (r + f)~ x = ((r + = 1 1 in H r 

where the average value T _1 is defined later. As shown 
in Ref.pj], the boundary condition of ^(r = 0) = holds 
when we choose E r = 0. This gives us the values of k as 
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A 2 cos 2 (0) + A 2 Z sin 2 ((/>). We can also get solu- 
tions with k being the negative of the expressions on 
the RHS. But since it is only the positive values which 
are compatible with the boundary condition that the 
states are localized near the surface, we focus our 
attention only on the positive solutions n±. To find the 
eigenvectors, we note that 

T Z 
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where in the last line we have used e l€ ^ = 
{A y coscj) + %A Z sin^)/^^. We then obtain the two ba- 
sis (normalised) eigenstates of H r given by \ip±) = 
A/0p(r, (j))\u±) where p(r, <p) = e~ K+r — e~ 
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and is a dependent normalisation factor. The con- 
vention chosen above is such that |^ + ) corresponds to 
the solution in which spin is pointing along the nega- 
tive r direction (ie outside the insulator) and the spin 
for points along the positive r direction (inside the 
insulator). Note that the direction of the real spin is de- 
termined by the angle (/>, which depends on the values of 
a, which determines the curve of the surface and also A y 
and A z which are material dependent parameters. It also 
varies as a function of <p , as we sweep the surface across 
the range of (j). But it is always parallel to the surface, 
as is expected from spin-momentum locking. 

Now, let us derive the effective 2D surface Hamiltonian 
in the l^zb) space. Any surface state \x) can be repre- 
sented as |x) = x+\^+) + X— IV 7 —)- We can now define 
a two component spinor x = (x+?X-) T an d define the 
effective Hamiltonian for x as 
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Note that + and — refers here to real spin-up and spin- 
down, but the quantisation axis, and hence the mean- 
ing of spin up and down, continuously changes along the 
hyperbolic surface. We shall essentially follow the same 
steps as in Ref.[19] to compute the effective Hamiltonian. 
To simplify the notation, we define #± = </>/2^F7r/4. We 
find that the diagonal elements (%j)±\H ( ] ) + H x \ip±) = 
and the off-diagonal elements are given by 
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As was done in Ref. [19], it is convenient to define a 
length variable, instead of the angle variable <p as 
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located just below the geometric surface. The limits of 
(j) —> tan -1 a,7r — tan -1 a correspond to I —> — oo, +00. 
Since the probability density must not change during this 
transformation, X±(l) is related to x±(4 > ) as 



(18) 



Hence, we can solve the eigenvalue equation H e ^x± — 
Ex± on the entire hyperbolic surface, by transforming 

h±x{4>) = e x {4>) to 



(19) 



4 



where fi± is now given by 

S)± = TAA t \di (X) ~ A xPx , (20) 

and A\ can be considered the effective velocity of the 
electron along the coordinate I which is effectively in- 
creased due to the effect of the Jacobian dependent sec- 
ond term in A^. We focus on eigenstates with energy 
E 



tions as 

X± = 



(A x k x ) 2 and obtain the two surface SOlU- 
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The factor of y/2 above ensures that the probablity den- 
sity current in the direction of I is =b^ as it should be 
for a free Dirac particle of momentum k and energy E. 
Hence, as expected, we get surface states satisfying the 
Dirac equation. However, unlike on a planar surface, the 
definition of Ah changes as we change I. Since these wave- 
functions are valid everywhere on the hyperbolic surface, 
it is clear that no backward scattering takes place any- 
where and that the transmission is unity along any path 
independent of the shape and size of the curvature of the 
surface. Thus, we generalise the earlier result [19] which 
concluded that there was no reflection at a 90° corner 
and find that there is no reflection even when the corner 
has any angle other than 90°. 



III. DENSITY OF STATES AND TUNNELING 
CURRENT FOR POLARISED STM 

A. Density of states 

Although there is no reflection on the curved surface, 
the density of states does change as a function of both 
the curvature and the sharpness of the corners. The local 
density of states (DOS) is defined as 

p s = f duj[f(uj)-f(uj^eV)]Y,^u-eV)\x ± \ 2 (22) 

and is clearly a function of the curvature parameter a and 
the sharpness parameter R through the wave-function 
X±. We will restrict ourselves to low biases and low 
temperatures and normalize the DOS by its value at 
(j) = tan -1 a. In Figs. (2a) and (2b), we show the DOS 
as a function of the angle parameter <fi which spans the 
surface. Note that the DOS shows a dip at <j) = ir/2 (or 
I = 0), which is the point of maximum curvature, both 
for fixed R in Fig. (2a) as well as for fixed a in Fig. (2b). 
Note also that in Fig. (2a), the range of <p depends on 
the curvature parameter a and increases as the curva- 
ture increases, whereas in Fig. (2b), the range is fixed 



7r/4 to 7r — tan 



3tt/4 for a = 1. 



from tan - < 

From both the curves, it is clear that in the limit of a 
planar surface (a — >• oo or R — » oo), the DOS is flat, as 
is expected for the usual planar TIs. 





FIG. 2. The DOS as a function of the angle parameter <j> along 
the surface, for fixed sharpness R = 1 and for three different 
curvatures, normalised to unity at cj) = tan -1 a, is shown in 
Fig. (2a), while in Fig. (2b), it is shown for three different 
values of the sharpness parameter and for fixed curvature a = 
1. 



B. Tunneling current for a polarized STM tip 

We now use the theory for the tunneling current of 
Dirac electrons on the surface of a TI, developed in 
Ref . [T8] . By using the Bardeen tunneling formula[2T], 
and explicitly computing the matrix element for overlap 
between the tip and the material wave- functions, they 
found that the tunneling current at zero temperature was 
given by 

I(Y) = Io\co\ 2 Pt[Pd + p z cos a + Pm sin a] 
x J dcuif(cu)-f(cu + eV)} 

with p d = Y J KEv-eV)\^ v \ 2 (l^nl) 
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elOlWlVcosCS-^) .(23) 
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Here Co is a constant which depends on details of the 
tip wave-function and p t is the density of states for the 
tip electron, assumed to be constant. However this was 
derived using a flat surface for the topological insulator 
with the spin quantisation axis of the electrons fixed to 
be in the ^-direction. In this basis, the tip wave-function 
was represented as 



■A 



cos(a/2) 
sin(a/2)e i(3 



(24) 



and the electron wave-function on the surface was given 

by 



1 



*h> • (25) 



But for the hyperbolic surface, the spin quantisation axis 
is continuously changing as a function of the curvature - 
i.e., as a function of the angle 0, due to spin- momentum 
locking. Hence, to express the tip wave- function in the 
basis of the electron spin quantisation axis, we need to 
rotate the spinor by <\> — tt/2 and obtain it in the 0-basis, 
the basis which is perpendicular to the surface. In this 
basis, the tip wave- function is given by 



Ccos(a/2) -iSsm(a/2)e^ 
-iScos(a/2) + Csin(a/2)e i(3 



(26) 



where C = cos(0/2 - 7r/4) and S = sin(0/2 - tt/4). Note 
that we have rotated the basis in the counter-clockwise 
direction about the x-axis, and have chosen the 0-axis 
to be pointing out of the TI - i.e., opposite the normal 
direction defined by r. 

Now to be able to use the expression for the current 
given in Eq(23j we need to rewrite the spinor given in 



Eq 26 in the form given in Eq 24 We find that ipf can 



be written as 



cos(a/2) 
sin(d/2)e^ 



(27) 



where 

cos a = — cos (j) sin a sin /3 + sin <\> cos a, 

tan $ = tan f3 sin <\> + cos <\> sec f3 cot a (28) 

and the overall (unimportant) phase of the spinor is given 

by 



tan£ 



tan(f )cos(/?)tan(§ - f) 
l + tan(f)sin(/3)tan(f - f) 



(29) 



The wave-function of the electron on the curved surface 
has been derived in Eqj2l] Comparing with the wave- 
function in Eq 25 1 we see that n v = —1 and r\ v = d=£. 



pute the tunneling conductance as 

G = Go\co\ 2 pt[pd + Pz cos a + p m sin a] 
with p d = ^2,S{E V - eV), 

v 

Pz=0 

and p m = - X J2 S(E„ - eV) cos0 T ■ (30) 



where Go = 2e 2 /h. p z = because k v = — 1 , which 
is a result of spin-momentum locking, due to which the 
spin of the electron lies along the surface perpendicular 
to the quantization axis. However, p& and p m are non- 
zero and both show non-trivial dependence on (or I). 
In Fig. (3), the tunneling conductance has been plotted as 
a function of I for various values of the curvature a with 
fixed sharpness parameter R (for fixed tip parameters 
a and /?). Note that unlike the case for a flat surface 
for which the STM conductance is constant, (reproduced 
here in the large a limit), here the STM conductance 
varies as it spans the curved surface and shows a dip 
precisely at I = 0. This dip depends on the curvature 
and increases as the curvature increases. 
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FIG. 3. Tunneling conductance G (normalised to unity at 
/ — > — oo) for fixed polar and azimuthal angles (a = and 
/3 — 0), as a function of /, the parameter along the surface, 
for R — 1 and for three different curvatures. 



Although pd in Eq 



30 



Now, we can use the expressions given in Eq 23 to corn- 



is independent of a and /?, p m is 
dependent on both and hence shows a non-trivial depen- 
dence on the polar and azimuthal angles of the STM tip. 
Hence, the tunneling conductance also has a non-trivial 
dependence on both a and /?, which has been shown in 
Figs. (4a) and (4b) respectively. As can be seen from the 
plots, the conductance shows a dip, not only at (j) = tt/2 
( or I = for a ^ 1), but also at a = 0, tt in Fig. (4a) and 
at p = tt/2, 3tt/2 in Fig. (4b). 



6 



(a) 




to cases where the convex side is filled with the TI. As 
explained in Ref.[19 , in this case, r < 0. But the Jaco- 
bian becomes ill-defined when \r\ > R due to the pres- 
ence of the term (— \r\ + Hence, we need to re- 
strict ourselves to \r\ < R. In that case, a very similar 
analysis works and the surface is described by the effec- 
tive Hamiltonian in Eq(20]with the difference that now 

A<j> = [A4 - ra 2 ((-|r| + (d^fi). The sign change 

compared with the earlier case implies that the velocity 
of the electrons is now reduced, and hence its amplitude 
and consequently, the DOS shows peaking at the corners 
instead of dips. This can also be measured by an STM 
tip. 

In conclusion, we have found the effective Hamiltonian 
on curved hyperbolic surfaces and shown that the Hamil- 
tonian can be written in terms of a continuous coordinate 
which varies along the surface. This clearly indicates that 
there is no back-scattering for any curvature. The sharp- 
ness of the edge can also be changed continuously with no 
backseat tering even at sharp edges. However, the DOS 
does change as we span the surface and there is a dip (or 
peaking) of the DOS at the edge for concave (convex) 
surfaces. We have also shown that the STM spectra of 
the Dirac electrons on the curved surface, as measured by 
a magnetized tip, shows unconventional and non-trivial 
dependence, not only on the parameter spanning the sur- 
face, but also on the polar and azimuthal angles of the 
tip. These measurements would provide clear evidence 
for the curvature of the surface. 



FIG. 4. In Fig. (4a), the tunneling conductance G (nor- 
malised to unity at <j) = tan -1 a) is plotted as a function 
of the polar angle a of the tip and 0, the angle that spans 
the surface, for fixed azimuthal angle ft — 0, fixed curvature 
a — 1 and fixed sharpness parameter R = 1. In Fig. (4b), 
the tunneling conductance G is plotted as a function of the 

azimuthal angle f3 and 0, for fixed polar angle a = 0.1 for the ACKNOWLEDGMENTS 

same curvature and sharpness. 



IV. DISCUSSION AND CONCLUSION 

Although we have restricted ourselves to concave hy- 
perbolae for the above calculations, they can be extended 



One of us (S.R.) would like to thank Sourin Das, 
Venkat Pai, Arijit Saha, Diptiman Sen, Vijay Shenoy and 
Abhiram Soori for useful discussions. We would partic- 
ularly like to thank Arijit Saha for bringing Ref.[l9] to 
our attention. 



[1] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005); ibid, Phys. Rev. Lett. 95, 146802 (2005); C. Wu, 

B. A. Bernevig, and S. C. Zhang, Phys. Rev. Lett. 96, 
106401 (2006); BA. Bernevig, T.L. Hughes, and S.C. 
Zhang, Science 314, 1757 (2006); J. E. Moore and L. 
Balents, Phys. Rev. B75, 121306 (2007); R. Roy, Phys. 
Rev. B79, 195321 (2009). 

[2] Konig, M., H. Buhmann, L.W. Molenkamp, T.L. Hughes, 

C. -X. Liu, X.L. Qi, and S.C. Zhang, J. Phys. Soc. Jpn. 
77, 031007(2008); Konig, M., S. Wiedmann, C. Brune, 



A. Roth, H. Buhmann, L. Molenkamp, X.-L. Qi, and S.- 
C. Zhang, Science 318, 766 (2007). 

[3] M. Buttiker, Science 325, 278 (2009); M. Z. Hasan and 
C. L. Kane, Rev.Mod.Phys. 82, 3045 (2010), X-L Qi and 
S-C Zhang, Rev. Mod. Phys. 83, 1057 (2011). 

[4] See, for instance, O. Heinonen, 'Composite fermions: a 
unified view of the quantum Hall regime', World Scien- 
tific, Singapore (1998). 

[5] C. -Y. Hou, E. -A. Kim, and C. Chamon, Phys. Rev. 
Lett. 102, 076602 (2009) ; J. Maciejko, C. Liu, Y. Oreg, 



7 



X. -L. Qi, C. Wu, and S. C. Zhang, Phys. Rev. Lett. 102, 

256803 (2009) ; A. Strom and H. Johannesson, Phys. 

Rev. Lett. 102, 096806 (2009) ; A. Bermudez, D. Patane, 

L. Amico, M. A. Martin-Delgado, Phys. Rev. Lett. 102, 

135702 (2009) ; J. E. Moore, Nature 464, 194 (2010); 

N. Goldman, I. Satija, P. Nikolic, A. Bermudez, M. A. 

Martin-Delgado, M. Lewenstein, I. B. Spielman, Phys. 

Rev. Lett. 105, 255302 (2010) ; R. Egger, A. Zazunov 

and A. Levy Yeyati, Phys. Rev. Lett. 105, 136403 (2010); 

Joseph Maciejko, Eun-Ah Kim and Xiao-Liang Qi, Phys. 

Rev. B 82, 195409 (2010); J. E. Vayrynen and T. Ojanen. 

Phys. Rev. Lett. 106, 076803 (2011); P. Virtanen and P. 

Recher, Phys. Rev. B 83, 115322 (2011). 
[6] S. Das and S. Rao, Phys. Rev. Lett. 106, 236403 (2011). 
[7] A. Soori, S. Das and S. Rao, Phys. Rev.B 86, 125312 

(2012). 

[8] T. Zhang, P. Cheng, X. Chen, J. Jia, X. Ma, K. He, L. 

Wang, H. Zhang, X. Dai, Z. Fang, X. Xie, and Q. Xue, 

Phys. Rev. Lett. 103, 266803 (2009). 
[9] P. Roushan, J. Seo, C. V. Parker, Y. S. Hor, D. Hsieh, 

D. Qian, A. Richardella, M. Z. Hasan, R. J. Cava and A. 

Yazdani, Nature 460, 1106 (2009). 



[10] J. Seo, P. Roushan, H. Beidennkopf, Y.S. Hor, R. J. Cava 

and A. Yazdani, Nature 466, 343 (2010). 
[11] K. Imura, Y. Takane and A. Tanaka, Phys. Rev. B84, 

195406 (2011). 

[12] K. Imura, Y. Yoshimura, Y. Takane and T. Kukui, cond- 

mat/1205.4878. 
[13] K. Imura and Y. Takane, cond-mat/1211.2088. 
[14] B. Zhou, H. Lu, R. Chu, S. Shen, and Q. Niu, Phys. Rev. 

Lett. 101, 246807 (2008). 
[15] Y. Jiang, F. Lu, F. Zhai, T. Low and J. Hu, Phys. Rev. 

B84, 165439 (2011). 
[16] R. Takahashi and S. Murukami, Phys. Rev. Lett. 107, 

166805 (2011). 

[17] D. Sen and O. Deb, Phys. Rev. B85, 245402 (2012); S. 

Modak, K. Sengupta and D. Sen, cond-mat/1203.4266. 
[18] K. Saha, S. Das, K. Sengupta and D. Sen, Phys. Rev. B 

84, 165439 (2011). 
[19] Y. Takane and K. Imura, J. Phys. Soc. Jpn, 81, 093705, 

(2012). 

[20] H. Zhang, C. Liu, X. Qi, X. Dai, Z. Fang and S. C. Zhang, 

Nature Physics 5, 438 (2009). 
[21] J. Bardeen, Phys. Rev. Lett. 6, 57 (1961). 



